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SUMMARY 

Three different computer codes developed in-house are described for application 
to hot composite structures. These codes include capabilities for: (1) laminate 

behavior (METCAN), (2) thermal/structural analysis of hot structures made from high 
temperature metal matrix composites (HITCAN), and (3) laminate tailoring (MMLT) . 
Results for select sample cases are described to demonstrate the versatility as well 
as the application of these codes to specific situations. The sample case results 
show that METCAN can be used to simulate cyclic life in high temperature metal matrix 
composites; HITCAN can be used to evaluate the structural performance of curved 
panels as well as respective sensitivities of various nonlinearities; and MMLT can be 
used to tailor the fabrication process in order to reduce residual stresses in the 
matrix upon cool-down. 


INTRODUCTION 

High temperature metal matrix composites (HT-MMCs ) are emerging as materials 
with potentially high payoffs in aerospace structural applications. Realization of 
these payoffs depends on the parallel and synergistic development of: (1) a tech- 

nology base for fabricating HT-MMC structural components, (2) experimental techniques 
for measuring their thermal and mechanical characteristics, and (3) computational 
methodologies for predicting their nonlinear behavior in complex service environ- 
ments. In fact, it might be argued that the development of computational methodolo- 
gies should precede the others because the structural integrity and durability of 
HT-MMCs can be computationally simulated, and the potential payoff for a specific 
application can be assessed, at least qualitatively. In this way, it is possible to 
minimize the costly and time consuming experimental effort that would otherwise be 
required in the absence of a predictive capability. 

Recent research at NASA Lewis is directed towards the development of a com- 
putational capability to predict the nonlinear behavior of HT-MMCs. This capability 
is in the form of stand-alone computer codes which are used to computationally 
simulate HT-MMC behavior in all its inherent scales. The simulation starts with 
constituents and the fabrication process and proceeds to determine the effects 
induced by the severe service loading environments. Three computer codes have been 
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developed to-date: (1) METCAN - Metal-Matrix Composite Analyzer, (2) HITCAN - High 

Temperature Composite Analyzer, and (3) MMLT - Metal Matrix Laminate Tailoring. 

These codes were developed in the order mentioned above. The amount of discussion 
allotted for each code reflects directly the experience we had with each code. Each 
of these codes simulates a specific area of HT-MMCs ♦ For example: (1) METCAN is 

structured to simulate HT-MMCs behavior at a point, (2) HITCAN is developed to 
analyze the behavior of hot structures with active cooling provisions (passages) made 
from HT-MMCs, and (3) MMLT is developed to concurrently tailor the constituent 
materials characteristics and the fabrication process for an a priori specified 
HT-MMC behavior such as minimum residual stresses upon cool-down. The primary 
objective of this article is to briefly describe these three computer codes and 
present illustrative results from their applications to simulate specific HT-MMC 
behavior. The secondary objective is to demonstrate what can be done with these 
three computer codes but not how it is done. Select references are cited for that 
purpose. 


METCAN - METAL MATRIX COMPOSITE ANALYZER 

The structure of METCAN, its simulation capabilities, and typical results to 
illustrate the applications of these capabilities, are summarized in this section. 

The structure of METCAN parallels the fabrication process of matrix composites. 
A typical fabrication process is schematically illustrated in figure 1. The simu- 
lation capability in METCAN is depicted schematically in figure 2. METCAN has the 
capability to predict all aspects of HT-MMC behavior, including the fabrication 
process by using only room temperature properties for the fiber and matrix. The 
formalism embedded in it, an initial version, and concept demonstration are described 
in reference 1. A detailed description of the micromechanics to repreresent the 
simulation at the constituent materials level is provided in reference 2. 

Fundamental to the computational simulation in METCAN is the introduction of 
the multifactor interaction model (MFIM ) to represent the various nonlinearities and 
their mutual interactions in the constituents. The equation form of the MFIM and 
reasons for its selection are summarized in figure 3. A discussion on its ability to 
represent constituent material behavior and the subsequent influence of this behavior 
on the response of structural components from HT-MMC is presented in reference 3. 

The use of METCAN to simulate metal matrix composite behavior from consitutent 
material properties at room temperature is summarized in table I. Comparisons of 
these properties with three-dimensional finite elements are shown in table II 
(ref. 4). METCAN simulation of the cyclic behavior of HT-MMCs is described in 
reference 4, where the influence of the interphase and limited comparisons with room 
temperature data are also described. Illustrative results from reference 4 are shown 
in figure 4. It can be observed in this figure that cyclic loading influences the 
MMC stress/strain behavior and degrades the stress at fracture. METCAN simulation of 
in-situ behavior and how this can be used to interpret composite-measured behavior 
are described in reference 5. Results from reference 5 on the effects of in-situ 
matrix strength on unidirectional MMC transverse tensile and shear strengths are 
shown in figures 5(a) and (b), respectively. Degradation of in-situ matrix strength, 
compared to its bulk state, substantially degrades matrix dominated strengths. The 
degradation is severest for the transverse tensile strength. The degradation for 
longitudinal tension is negligible but it is considerable for longitudinal compres- 
sion. Therefore, transverse tensile tests will be the most sensitive to experi- 
mentally detect any degradation. 
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Corresponding results for the development of an interphase between fiber and 
matrix, or weakening of the interfacial bond, are shown in figures 6(a) and (b) for 
transverse tensile and shear strengths, respectively. The interphase properties or 
the interfacial bond strength were assumed to be ratios of the matrix undegraded 
strength, while the in-situ matrix was assumed to retain its bulk state properties. 
It can be seen in figure 6 that the interphase degrades the stresses at fracture but 
does not degrade the stress/strain curve. There are other conditions of the in-situ 
matrix which can influence matrix strength as well. Some of these are shown in 
figure 7 (ref. 4). The annealed yield with no interfacial bond corresponds to very 
low values measured for these composites. In this context, METCAN can also be used 
to interpret measured data. 

METCAN simulation of life under thermal cycles is shown in figure 8, where 
available measured data (ref. 6) is also shown for comparison. This very good 
comparison is considered adequate to demonstrate the simulation capability of METCAN 
for cyclic life. 


HITCAN - HIGH TEMPERATURE COMPOSITE ANALYZER 

HITCAN combines METCAN with a noncommercial finite element code, MHOST, and a 
dedicated mesh generator. The code is stand-alone and stream-lined for the thermal/ 
structural analysis of hot metal matrix composite structures. A schematic of the 
code's structure iB shown in figure 9, with its capabilities summarized in table III. 
An extensive description of HITCAN, including a variety of sample cases to illustrate 
its computational capabilities, are found in reference 7. The results for a curved 
panel are included herein as specific examples. 

The panel geometry, laminate lay-up, and loading conditions (thermal, mech- 
anical) are shown in figure 10. The buckling evaluation results are shown in fig- 
ure 11. The buckling load decreases with fiber degradation (fiber diffusion into the 
matrix) and with temperature (material degrades due to temperature). The decrease, 
due to temperature, is substantial (about 30 percent). The panel vibration fre- 
quencies, displacement, ply and constituent stresses are summarized in figure 12. 

This figure illustrates the versatility (breadth and depth) of HITCAN to simulate MMC 
structural behavior and evaluate their adequacy in hot structures applications. 

HITCAN can also be used to perform sensitivity analyses. Results from sensitivity 
analyses for the curved panel are shown in figure 13. 

Collectively, the results summarized demonstrate that the complex behavior of 
MMC structures can be simulated at all its scales by using an integrated computer 
code such as HITCAN. 


MMLT - METAL MATRIX LAMINATE TAILORING 

The Metal Matrix Laminate Tailoring computer code consists of METCAN with a 
suitable optimizer. This code has been recently initiated and only preliminary 
results have been obtained to date. The code is described in detail in reference 8. 
Herein, representative results are included to illustrate its application. The 
results are for a graphite fiber/copper-matrix composite where the fabrication 
process (temperature and consolidation pressure) histories are tailored to minimize 
the residual stresses in the matrix during cool-down from consolidation to room 
temperature. Constraints (conditions) were imposed in the tailoring procedure that 


3 



the residual stresses will not exceed the corresponding matrix strength during 
cool-down. 

The tailored fabrication process is shown in figure 14. The microstress 
developed in the different matrix regions (A, B, or C, fig. 3) are shown in fig- 
ure 15. These results demonstrate that MMLT can be used to process HT-MMCs for 
desired matrix stress magnitudes during the cool-down process as well as the 
sensitivities of the various parameters that influence the optimum fabrication 
process. For example, in the present study, the consolidation pressure history 
appeared to be one of the more important parameters that influence the optimum 
fabrication process. 


CONCLUSIONS 

Three computer codes have been developed and are described which can be used to 
simulate the complex behavior of hot structures made from high temperature metal 
matrix composites (HT-MMCs). These codes are: (1) Metal Matrix Composite Analyzer 

(METCAN), (2) High-Temperature Composite Analyzer (HITCAN) , and (3) Metal Matrix 
Laminate Tailoring (MMLT). Results from each code for select sample cases are 
included to illustrate the capabilities of each code. The results from METCAN are 
for laminate properties, in-situ strength, interphase effects, and cyclic load 
effects. The results from HITCAN are for a curved panel subjected to thermal and 
mechanical loads and including various nonlinearities. HITCAN results for sensit- 
ivity analyses are also included. The results for MMLT are for tailoring the 
fabrication process in order to minimize the residual stresses in the matrix during 
cool-down. Collectively, the results from these sample cases demonstrate that 
computer codes can be developed to computationally simulate composite hot structure 
at all its scale levels. 
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TABLE I. - NOMINAL CONSTITUENT PROPERTIES AND PREDICTED COMPOSITE 

PROPERTIES AT ROOM TEMPERATURE 


(a) SiC fiber 


Fiber density, p f , lb/in, 3 0.11 

Fiber elastic modulus, E f , Mpsi . 62 

Fiber Poisson's ratio, U f , in. /in 0.3 

Fiber shear modulus, G f , Mpsi 23.8 

Fiber coefficient of thermal expansion, 0f f , ppm 1.8 

Fiber melting temperature, T Mf , °F 4870 

Fiber tension strength in direction 11, s fUT > ksi 500 

Fiber composite strength in direction 11, S fllc , ksi 650 

Fiber tension strength in direction 22, S f22T , ksi 500 

Fiber composite strength in direction 22, s t22c r ksi 650 

Fiber shear strength in direction 12, S ksi 300 

Fiber diameter, D f , mils 5.6 


(b) Ti-15-3-3-3 Matrix 


Matrix density, p , lb/in. 3 0.172 

Matrix elastic modulus, E^, Mpsi 12.3 

Matrix Poisson's ratio, V , in. /in 0.32 

m . _ 

Matrix shear modulus, G^, Mpsi ..... 4.7 

Matrix coefficient of thermal expansion, O^, ppm . . ^ - 4.5 

Matrix melting temperture, T , °F 1800 

Matrix tension strength, S , ksi 130 

Matrix composite strength, S^ r ksi 13° 

Matrix shear strength, S , ksi 91 


(c) Composite 


Composite density, p^ f lb/in. 3 * 0.147 

Composite eleastic modulus in direction 11, E , Mpsi . 31.2 

Composite elastic modulus in direction 22, E c22 , Mpsi 18.5 

Composite elastic modulus in direction 33, E c33 , Mpsi 18.5 

Composite shear modulus in direction 12, G cl2 , Mpsi 7.7 

Composite shear modulus in direction 23, G c23 , ^ * 9 

Composite shear modulus in direction 13, G , Mpsi 7.7 

Composite Poisson's ratio in direction 12, V , in. /in 0.31 

Composite Poisson's ratio in direction 23, V c23 r in. /in 0.33 

Composite Poisson's ratio in direction 13, V cl3 , in. /in 0.31 

Composite coefficient of thermal expansion in direction 11, 

a cll , ppm 2 * 9 

Composite coefficeint of thermal expansion in direction 22, 

<* c22 ' 3-2 

Composite coefficient of thermal expansion in direction 33, 

Q c33' 3 ' 2 
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TABLE II. - GRAPHITE P-IOO/COPPER ROOM TEMPERTURE MECHANICAL AND THERMAL PROPERTIES: COMPARISON OF METCAN 

PREDICTIONS AND THREE-DIMENSIONAL FINITE ELEMENT ANALYSIS 
Property type Fiber volume ratio 

’ 0.063 I 0.2234 I 0.466 I 0 . 622 



METCAN 

3-DFEM 

METCAN 

3-DFEM 

METCAN 

3-DFEM 

METCAN 

3-DFEM 

E |n' mp0i 

23.7 


37.3 

36.7 

58.4 

58.0 

72 . 0 

71.8 

E | 22 ' m V Bi 

13.8 


10.2 

11.0 

6.5 

6.9 

4.6 

5.0 

E |33' mpSi 

13.8 


10.2 

11.0 

6.5 

6.9 

4.6 

5.0 

°fl2' mpSl 

5.8 

6.0 

4.5 

4.8 

3.2 

3.4 

2.5 

2.6 

°*23' mpBl 

5.6 

5.8 

4.2 

4.5 

2.8 

2.7 

2.1 

1.8 

G| 13 , mpsi 

5.8 

5.8 

4.5 

4.5 

3.2 

2.7 

2.5 

1.0 

V tl2' ln ' /ln * 

0.29 

0.30 

0.28 

0.30 

0.25 

0.29 

0.24 

0.28 

V t2 y in./in. 

0.30 

0.30 

0.30 

0.26 

0.30 

0.22 

0.30 

0.19 

in. /in. 

0.30 

0.30 

0.30 

0.26 

0.30 

0.22 

0.24 

0.22 










a fll' ppm 

6.8 

6.6 

3.3 

3.7 

1.0 

1.40 

0.18 

0.42 

“|22' PPm 

10.4 

10.5 

11.0 

10.9 

11.2 

11.1 

11.1 

10.9 

a |33' PE>m 

10.4 

10.5 

11.0 

10.9 

11.2 

11.1 

11.1 

10.9 

K * , Btu in./°F hr in. 2 

19.7 

20.6 

20,6 

21.6 

22.0 

23.2 

22.8 

24.0 

* 2 
Btu in./°F hr In. 

15.6 

17.3 

11.8 

13.5 

7.8 

8.8 

5.8 

6.3 

1 22' , 









K* , Btu in./°F hr in. 

15.6 

17.3 

11.8 

13.5 

7.8 

0.8 

5.8 

6.3 

1 ^ J 




























TABLE III. - HITCAN CAPABILITIES FOR COMPOSITE MATERIALS 


Type of analysis 


Type of structure 


Beam 

Plate 

Ring 

Curved 

panel 

Builtup 

structure 

Static 

Tested 

Tested 

Tested 

Tested 

Tested 

Buckling* 






Load stepping 






Modal (natural vibration modes) 






Time domain 






Loading: 





Tested 

Mechanical 

Tested 

Tested 

Tested 

Tested 

Thermal 

Tested 

Tested 

Tested 

Tested 

Tested 

Cyclic 






Impact 






Constitutive models 0 : 






P = constant 

Tested 

Tested 

Tested 

Tested 

Tested 

P * f (T) (temperture dependence) 

P - f (0) (stress dependence) 

P = f{0) (stress rate dependence) 






P = f(t) (creep) 





Tested 

P = f(T ,0,0) (combination) 

Tested 

Tested 

Tested 

Tested 

P » f (T, 0,0 f t) (creep combination) 






Fiber Degradation 

Tested 

Tested 

Tested 

Tested 

Tested 

Fabrication induced stresses 

Tested 

Tested 

Tested 

Tested 

Tested 

Ply orientations' 1 (arbitrary) 

Tested 

Tested 

Tested 

Tested 

Tested 


“Tested one buckling mode. 
b Tested four vibration modes. 

Constitutive model notation: P, material properties; T , temperature; 

O, stress; <7, stress rate; t , time. 

d Tested three ply orientations: unsymmetric (0/+45/90); symmetric (0/45 ) s ; 

balanced (0/90 ) ^ 
























Mix 


Heat 


Roll 


Matrix 

powder 

+ 

Teflon 

binder 

+ 

Stoddard 

solution 



Lay up and 
hot press 


Powder 

cloth 



Composite 


Figure 1 .-Metal-matrix composite fabrication process. 
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"Synthesis" 




Global 

structural 

analysis 



Finite element 


\ 

Laminate 



Composite 

micromechanics 

theory 



Composite 

micromechanics 

theory 


Constituents 

Figure 2— Integrated approach to metal- matrix composite analysis. 
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Matrix 



Notations: 

P - property; T - temperature; S - strength; R - metallurgical reaction; N - number of cycles: 
t - time; over dot - rate; subscripts: O - reference; F - final; M - mechanical; T — thermal 

Figure 3.— Assumed multifactor interaction relationship to represent the various factors which influence in situ constituent 
materials behavior 



0 .24 .48 .72 .96 1.20 0 .24 .48 .72 .96 1.20 

Strain EPS XX, percent 

Figure 4. — Metcan predictions for thermo mechanical cycling of MMC (SiC/Ti 15 333 at .36 FVR). 
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In-situ matrix strength 
A: 100% of reference value 
B: 50% of reference value 
C: 25% of reference value 




Figure 5.— Influence of in-situ matrix strength on the response of SiC/Ti- 15-3-3-3 unidirectional composite. 


Interphase properties 
A: Same as matrix 
B: 50% of matrix 




Figure 6.— Influence of interphase property degradation on the response of SiC/Ti- 15-3-3-3 unidirectional 
composite. 



Figure 7 —Unidirectional graphite/copper MMC transverse 
strength bounds. 
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Number of thermal cycles 

Figure 8 — Metcan simulates thermal fatigue 
strength of unidirectional SiC/Ti composite. 
(FVR .35; T = 70 °F). 






Fixed-free curved panel under bending and uniform temperature 
loadings for (SiC/Ti- 15-3-3-3, 0/±45/90); 

0.4 fiber volume ratio 

Geometry 



Loading 



Load steps 


Figure 10 —Geometry and loading for curved panel. 


Fixed-free curved panel under external pressure and uniform 
temperature loadings for (SiC/Ti-15-3-3-3, 0/±45/90); 

0.4 fiber volume ratio 

Geometry and loading 



(1) Under mechanical loading only = 38053 psi 

(2) With fiber degradation, under mechanical loading only 
- 35437 psi 

(3) Under thermo-mechanical loading = 25000 psi 

Figure 1 1 .—Buckling analysis for curved panel under thermo- 
mechanical loading. 
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Fixed-free curved panel under bending and uniform temperature loadings for (SiC/TT-1 S-3-3-3, 0/±45/90); 0.4 fiber volume ratio 


Geometry 


Loading 



1 2 3 

Load steps 



Displacements, 
curved edge center 



IX. 

0 

© 



_ p iy 

Inlerphase (B) 

Matrix^ A) 


1 2 3 

Load steps 


Stresses, curved edge center, ply 1 
0 r 

P -25 

-Ply 


Interphase (B) 
interphase (C) 



-** x r Matrix (C) 
' Lc- 


1 2 3 

Load steps 

Figure 12— Base case results for curved panel. 
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I 1 I 


1 2 3 

Load steps 
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Fixed-free curved panel under bending and uniform temperature loadings for (SiC/Ti-15-3-3-3, 0/±45/90); 0.4 fiber volume ratio 

Response: at load step 3 


Geometry and loading 
' '» 2000 


Fabrication-induced stresses 


Location of 

displacement 

response 




<D 

S 1500 

! ,000 
m 

E 500 
B 

i* 



1 2 3 

Load steps 


Fabrication- 

induced 

stresses 

Disp. 

[free 

end 

center] 

(inch) 

No 

-0.0157 


-0.0130 


Stresses, ply 1 
[fixed end center] i 

(ksi) 

Longitudinal 

Transverse 

Shear 

40.9 

-48.6 



Fiber degradation 


Degradation 


Stresses, ply 1 
[fixed end center] (ksi) 



center] 

(inch) 

Longitudinal 

Transverse 

Shear 

No 

-0.0157 

40.9 

-48.6 

1.4 

Yes 

-0.0170 

26.3 

436.9 

1.1 


Ply orientations 


Orientation end 
center] 
(inch) 

(0/±45/90) -0.0157 

(0/45) s -0.0115 

(0/90) 8 -0.0115 



Stresses, ply 1 
[fixed end center] i 

[ksi) 

Longitudinal 

Transverse 

Shear 

40.9 

-48.6 

1.4 

20.6 

-48.4 

0.3" 

18.4 

-49.0 

0.0 


Constitutive relationships 
(nonlinear multi-factor interaction model) 


Relationship 

Disp. 

[free 

end 

center] 

(inch) 

P = constant 


P = f(T) 

Temp, dependence 

-0.0153 

P = f(o) 

Stress dependence 

-0.0138 

P « f(*) 

Stress rate dependence 

-0.0135 

P = f(T. o, 6) 
Combination 

-0.0157 


Notation: 

P = Material property 
a = Stress 


Stresses, ply 1 
[fixed end center] 

(ksi) 

Longitudinal 

Transverse 

Shear 

364 

-57.8 

1.1 

44.3 

-54.9 

1.3 

34.6 

-50.2 

1.1 

36.0 

-57.8 

1.1 

40.9 

-48.6 

1.4 


T = Temperature 
6 = Stress rate 


Figure 13. — Sensitivity analysis. 


15 









































16x10 z 



Figure 14. — Optimum and current cool-down phases for PiOO/copper. 


Longitudinal stress Transverse stress 



Figure 15. — Matrix microstresses developed during the cool-down phase of PiOO/copper. 
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